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SUMKAHY 

4 

A  ■imple  computational  method,  baaed  on 
the  simplex  algorithm  of  linear  programming , 
it  proposed  for  the  following  problem: 

"Consider  a  network  (e.g.,  rail,  road, 
oosaaunicatlon  network)  connecting  two  given 
polnta  by  way  of  a  number  of  intermediate 
points,  where  each  link  of  the  network  has  a 
number  assigned  to  It  representing  its  capacity. 
Assuming  a  steady  state  condition,  find  a  maxi— 

j 

»al  flow  from  one  given  point  to  the  other. 
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COMPUTATION  OP  MAXIMAL  FLOWS  LN  NETWORKS 

by 

D.  R.  Fulkerson 
0.  B.  Dant^ig 

'Hie  linear  programming  formulation  of  the  network-flow 
problem  given  In  fl]  Is  not  only  useful  a3  a  theoretical  tool : 
When  suitably  Interpreted,  It  provides  In  addition  a  simple 
and  efficient  hand— computing  scheme.  It  Is  our  aim  In  section  I 
to  describe  this  computation  for  networks  with  capacity  con¬ 
straints  on  area  only.  It  Is  unnecessary  to  use  programming 
terminology  In  this  description,  and  so  we  dispense  with  It. 

A  direct  proof  will  be  given  In  section  II  that  the  computational 
method  assures  one  of  finding  a  maximal  flow  In  an  arbitrary 
ne  twork . 

I .  Description  of  Method 

The  problem  (see  Llj,  [?] )  may  be  stated  as  follows: 

One  Is  given  a  connected  network  of  arcs  and  nodes  with  two 
distinguished  nodes,  called  source  and  sink  respectively . 

All  other  nodes  are  called  Intermediate.  Each  arc  In  the 
network  has  associated  with  it  a  positive  Integer,  Its  flow 
capacity.  The  direction  of  flow  Is  assumed  to  be  away  f rom 
the  source  In  arcs  having  the  source  as  endpoint,  and  Into 
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the  sink  in  sink  arcs;  otr.erwlse  ;io  flow  direction  is  specific  1. 
Subject  to  the  conditions  that  the  flow  in  an  arc  does  not 
exceed  it9  capacity,  and  that  the  total  flow  Into  any  inter¬ 
mediate  node  is  equal  to  *  he  flow  out  of  it,  It  13  desired 
to  find  a  maximal  flow  in  the  network,  i.e.,  a  flow  which 
maximizes  the  sum  of  the  flows  in  source  (or  oink)  arcs. 

In  order  to  illustrate  the  computing  method,  let  us  consider 
a  simple  example.  Suppose  we  have  the  network  of  Pig.  1  with 
source  A, 


sink  B,  arid  arc  capacities  as  indicated.  To  start  out,  select 
two  trees*  of  arcs,  one  branching  out  from  A,  the  other  from  D, 


A  tree  is  a  connected  linear  graph  without  circuits;  that  is, 

U  ere  is  one  and  only  or.e  chain  of  arcs  Joining  each  pair  of  nodes. 
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ao  that  every  intermediate  node  In  reached  i  y  Ju3l  one  of  the 
trees.  Call  these  TA  and  Tg.  For  example,  Tg  might  consist 
of  B  alone,  and  right  bo  the  net  AC,  CD,  DE ,  E F  (see  Pig. 
la).  Notice  that  since  the  network  13  connected,  It  Is  always 
possible  to  select  two  3uch  trees.  Next  Introduce  ary  arc 
which  leads  from  to  Tg .  There  will  then  le  Just  one  chain 
from  A  to  B;  flow  as  mucr.  as  possible  along  this  chain.  In 
the  example,  EB  Is  such  an  arc,  and  we  have  then  the  flow 
diagram  of  Fig.  2  with  all  of  the  arcs  AC, 


(P) 

Fig  2 


CD,  DE  saturated.  Select  any  one  of  these  saturated  arc3, 

say  AC,  and  place  some  identifying  mark  or,  It  for  future  reference. 

In  Fig.  2  we  have  used  a  bar  (|  }  ;  this  symbol  will  he  used 


It  Is  clear  that  arcs  may  be  removed  until  a  tree  Is  left. 
There  Is  then  a  unique  chain  Joining  A  and  B.  Elimination  of  any 
arc  of  this  chain  gives  two  trees  of  the  kind  described. 
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throughout  to  designate  a  subset  of  the  saturated  arcs.  Now 

observe  that  If  the  barred  arc  is  dropped  from  the  picture, 

we  again  have  two  trees  (Pig.  2a )  -  A  and  Tp  -  EB,  ED,  DC,  EP 

Again  introduce  any  pnbarred  arc  leading  from  to  Tp,  say 

AP.  This  creates  a  flow  along  the  chain  AF,  PE,  EB  of  10  units 

and  saturates  each  arc  of  this  chain.  Select  one  of  these, 

say  AP,  and  bar  it.  We  now  have  the  following  diagram  (Pig.  3), 


Fig  3 


and  we  have  achieved  a  flow  of  20.  Dropping  barreu  arcs  gives 
the  same  trees  as  in  Pig.  2a.  Introduce  arc  AE  fi-om  to  T  . 
This  leads  to  a  situation  we  have  sot  met  previously  in  that 
the  chain  thus  constructed,  namely  AE,  EB,  car.  not  ta.<e  any 
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more  f  low  because  EP,  thou./!  unbarred,  Is  it  its  ot  ,  •  r  1  !  ••  1  t . 
Liar  En  an  i  leave  AE  l;  with  a  flow  of  sere,  '  ot  itnl/.r  F'V.  4 
with  new  trees  33  si  own  Jn  Fig.  4a. 


Introduce  DB  to  get  the  chain  AE,  ED,  DB.  This  time  we  can 
get  an  Increase  even  though  DE  Is  saturated,  since  the  flow 
in  Pig.  4  Is  from  D  to  E.  Thus  If  the  flow  from  A  to  E  Is 
Increased  by  t  >  0,  the  flow  from  D  to  E  must  be  decreased 
by  A  ,  and  the  flow  from  D  to  B  Increased  by  -  (see  Pig.  4b) 
In  order  to  preserve  the  conservation  equations  at  E  and  D. 
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The  largeet  pooalble  value  of  £  la  10,  since  the  capacity  of 
DB  (and  of  AE)  la  10.  This  cancels  the  flow  from  D  to  E.  Bar 
DB  and  proceed. 

A  repetition  of  atepa  of  the  kind  described  produces  the 
sequence  of  flows  depicted  In  Flga.  3-8  below. 


Fig  5 


Fig  6 
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and  there  are  no  more  arcs  to  Introduce  from  TA  to  TQ.  At 
this  stage  examine  the  barred  arcs  connecting  nodes  of  TA  to 
those  of  Tg.  If  the  flow  in  each  of  these  is  In  the  "right" 
direction,  that  is,  from  to  T0,  an  optimum  has  been  reached, 
because  these  arcs  now  constitute,  In  the  original  network,  a  cut 
separating  A  and  B  (i.e.,  every  chain  Joining  A  and  B  contains 
one  of  these  arcs),  and  the  total  flow  through  the  network  is 
equal  to  the  sum  of  the  capacities  of  the  arcs  forming  the  cut. 

Hence,  since  it  is  clear  that  no  flow  can  exceed  the  sum  of  capa¬ 
cities  of  arcs  forming  a  cut,  the  flow  is  maximal  in  this  case. 

If,  on  the  other  hand,  the  flow  in  one  of  the  barred  arcs  which 

Join  Ta  to  Tg  Is  in  the  wrong  direction,  an  increase  in  total 

flow  may  possibly  be  obtained  by  decreasing  the  flow  in  this 

arc.  To  see  this,  notice  that  the  arc  in  question,  together  with  arcs 
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of  ar.'i  T^,  will  form  a  (unique)  chain  Joining  A  mj  E 
which  night  look,  for  example,  like  this 

©-T7-0 -rrOtTTOrrO-Tr® 


the  arrowheads  denoting  directions  of  flow.  Now  we  nay  subtract 
i  >  0  frorri  the  flow  In  all  arcs  having  the  sane  direction  as 
the  barred  arc,  add  €  to  the  flow  In  the  other  arcs.  Since 
these  latter  arcs  will  of  necessity  Include  the  source  (and  sink) 
arc  In  the  chain,  the  total  flow  will  be  Increased  provided 
there  Is  any  slack  In  the  chain,  that  Is,  If  none  of  the  arcs 
where  appears  Is  already  saturated  (though  unbarred  In  our 
procedure).  In  any  case,  locate  an  arc  In  the  chain  which 
determines  the  maximum  value  of  t  and  place  a  bar  on  It.  This 
either  gives  rise  to  two  new  trees  or,  If  the  rosltlor  of  the 
bar  has  rot  h^cn  changed,  reverses  the  flow  completely  In  fhe 
original  arc . 

In  the  example  (see  Fig.  k )  the  arc 3  DP,  DP ,  AE,  AP,  CF 
Joining  nodes  of  to  those  of  all  have  flows  Ir,  the  right 
direction.  Hence  the  flow  shown  in  Pig.  8  Is  maximal,  and  con¬ 
sequently  these  arcs  form  a  minimal  cut,  l.e.,  a  out  the  sum  of 
whose  arc  capacities  is  minimal. 
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It  is  clear  from  the  preceding  discussion  that  the 
prop>o3ed  computation  provides  a  nondecreasing  sequence  of  flows. 
All  practical  experience  with  the  method  indicates  that  a  maximal 
flow  is  reached  very  quickly.  There  i3,  however,  the  theoretical 
possibility  that  an  Infinite  sequence  of  flows  having  the  same 
value  may  be  obtained.  In  the  rest  of  the  paper  we  shall  indicate 
one  way  of  removing  this  difficulty.  The  idea,  not  a  new  one, 
is  to  perturb  the  capacities  of  the  area  Blightly  30  that  the 
decision  at  each  stage  as  to  the  placement  of  the  bar  becomes 
unique,  thus  assuring  a  finite  process.* 


II .  Proof  of  Convergence 

Given  an  arbitrary  network  N,  let  a^,...,a  be  a  listing 
of  the  arcs  in  some  order  and  let  ci»*'*>cn  be  their  capacities. 

At  each  stage  of  the  computation,  the  arcs  a^,  1-1,..., n; 
are  divided  into  four  mutually  exclusive  classes  T^,  Tg,  S,  Z 
where 


» 

From  the  linear  programming  point  of  view  the  method  thus 
far  corresponds  to  the  original  simplex  algorithm  without  modi¬ 
fications  to  account  for  degeneracy,  with  this  difference.  When 
"ties"  occur  in  the  method  [l.c.,  when  a  bar  may  be  placed  in 
one  of  several  positions  at  some  step),  we  have  not  Indicated 
a  unique  choico,  but  rather  have  used  language  like  "pick  any 
one..."  If  we  interpret  this  to  mean  "choose  at  random  among 
the  poasibil lties , "  then  we  are  at  least  assured  that  with 
probability  one  the  computation  terminates,  as  we  shall  see. 
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lo  a  tree  branching  out  from  A, 

7^  1 3  a  tree  branching  out  from  3, 

S  consists  of  the  barred  arcs, 

7.  contains  all  other  arcs, 

and  each  node  is  an  endpoint  of  some  arc  of  or  Tg.  Call 
such  a  division  a  basic  partition  of  N. 


Theorem :  For  a  given  basic  partition  of  N,  let  arbitrary  flows 

be  assigned  to  arcs  of  S  and  Z.  Then  the  conservation  equations 

uniquely  determine  flows  in  the  arcs  of  TA  and  Tg.  Moreover, 

if  xi>  • ‘ are  the  amounts  of  the  assigned  flows,  then  the 

m 

amount  of  flow  in  each  arc  of  N  is  equal  to  V  6  x,,  where 

j-i  J  j 


°j  -  1,  0,  or  -1. 

The  proof  can  be  made  clear  by  considering  the  following 
example: 


Fig  9 
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wnere  the  clotted  flows  x^,...,Xg  are  assigned. 

Orient  each  arc  of  the  two  trees  In  some  fasr.lon,  say 
away  from  A  In  and  away  from  E  In  Tg  (as  shown  In  Pig.  9), 
the  convention  being  that  a  positive  number  on  an  arc  represents 
flow  In  the  same  direction  as  lt3  orientation.  Thus  If 
x,  y,  ?  are  the  flows  In  arcs  AC,  CE,  CF,  respectively ,  the 
equation  at  C  Is  x  —  y  —  z  -  0. 

That  the  flows  In  the  arcs  of  each  tree  are  uniquely 
determined  can  be  seen  by  solving  the  equations  recursively  in 
the  following  manner.  Start  with  an  outermost  node  In  one  of 
the  tree3,  E  for  example.  There  is  then  only  one  unknown  flow 
Into  E,  ne.oce  lt3  value  can  be  found;  here  it  Is  x^  —  x^. 
Similarly  determine  the  flows  in  all  arcs  branching  out  of  F , 

0,  H,  and  then  work  backward  in  the  tree. 

For  the  last  assertion  of  the  theorem,  it  Is  convenient 
to  solve  tne  equations  in  another  way.  Select  any  one  of  the 
assigned  flows  x^.  There  are  two  cases  to  consider.  Either 
the  arc  with  flow  Joins  two  nodes  of  the  same  tree  or  not. 
Suppose  It  Joins  two  nodes  of  T^,  as  EF  does,  for  Instance. 

Since  there  Is  a  unique  chain  Joining  E  and  P  1*~>  ,  rarely 

FC,  CE,  the  flow  xr  In  EF  can  he  ta/er  nare  of  by  assigning 
♦  x -  appronriate ly  to  the  arcs  of  tf  Is  cnaln.  Here  the  assign¬ 
ment  1 8  ♦  x_  to  CF,  —  x.  to  CE .  O,  the  other-  ^a^q.  If  the  arc 
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in  question  Joins  TA  to  T^,  as  does,  consider  the  chain 
Joining  A  and  B  and  having  this  arc  as  one  of  its  links.  Here 
it  is  AC,  CE,  El,  IL,  LB,  and  the  assignments  are  x^  ,  x,,  x^ ,  — . x^, 
— x^ ,  respectively.  If  we  continue  in  this  fashion,  it  is  evident 
that  a  solution  to  the  system  of  equations  will  Le  built  up 
having  the  form  described  in  the  theorem.  For  the  example, 
the  solution  Is  Shown  In  Pig.  10. 


Fig  io 


In  the  algorithm  outlined  In  section  I,  arcs  of  S  are 
saturated  and  arcs  of  Z  have  zero  flow.  Hence  we  may  state 
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Co rollary : 
given  step. 


Suppose  S  consists  of  arcs  a.  , . . .  ,a.  at  any 

X1  xs 

Then  the  flow  In  each  arc  of  N  Is  given  by  a 

B 


sum  of  the  form  T  6  c.  ,  where  6.  *1,0,  or  —1. 

k-1  xk  xk  xk 

Let  us  examine  In  more  detail  how  ties  occur  In  the 
decision  as  to  which  arc  to  put  into  S  In  going  from  one  step 
to  the  next  in  the  computation.  Each  step  is  of  one  of  the 
following  two  types: 

(a)  a  flow  from  TA  to  Tg  is  Introduced  along  some  arc 
a  of  Z,  thus  shifting  one  of  the  arcs  of  or  Tg,  or  a  Itself, 
into  S; 


(b)  a  flow  from  Tg  to  TA  along  some  arc  b  of  S  is  decreased 
again  shifting  one  of  the  arcs  of  or  Tg  into  S,  or  leaving 

b  in  S  with  its  direction  of  flow  reversed. 

It  suffices  to  consider  (a),  since  (b)  Is  essentially  no 
different.  We  have  then  a  chain  from  A  to  B 


consisting  of  a  and  arcs  of  TA  and  Tg,  with  flows  In  one 
direction  or  the  other  in  each  arc  of  the  chain,  the  direction 
of  flow  in  the  source  (sink)  arc  being  away  from  A  (into  B). 


Suppose  the  arcs  In  the  chain  are  a 

and  let  S  consist  of  a.  ,  . . . ,  a. 

*1  x! 


1*1' 
J1  J2 


•  aj 


* ». 


m 
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By  the  corollary,  the  magnitude  of  flow  In  a  la  given  bv 

‘  l 

3  1 

5,JJ  ;  .  It  follows  that  we  may  Increase  the  flow  In  arc 

k-1  Ak  1k 

a  by  an  amount 

B  . 

<r  m  m  1  n  ( c  ,  +  7  6  c .  ) 

J,  b  b 

where  the  plus  sign  Is  taken  if  a.  has  flow  directed  toward 

Ji 

A,  the  minus  sign  otherwise.  Since  the  c,  are  integers,  it 

h 

is  clear  that  if,  Instead  of  using  the  original  capacities, 
we  had  altered  them  to  begin  with  by  setting 

C1  “  C1  +  10_1  '  1  ■  1 , • • • ,n  , 

the  minimum  would  be  achieved  for  Just  one  TVien  at  no  3tage 

do  ties  occur,  and  consequently  arcs  not  in  S  will,  at  each  step, 
be  unsaturated.  This  means  that  the  total  flow  through  N  is 
Increased  with  each  Iteration,  and  hence,  by  the  theorem,  no 
basic  partition  having  the  same  directions  of  flow  In  arcs  of 
3  car  reoccur.  Since  there  are  only  a  finite  number  of  basic 
partitions,  the  procedure  terminates  for  the  perturbed  problem. 
Thus  we  must  arrive  at  a  flow  In  which  all  area  of  S  leading  from 
Ta  to  Tp  form  a  cut  and  are  In  the  right  direction;  hence  the 
flow  1 8  maximal.  Once  an  optimum  has  been  attained,  a  solution 
to  the  original  problem  Is  given  by  rounding  each  arc  flow  to 
the  nearest  Integer. 
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Notlce  that  the  theorem  (see  [l],  [?J  )  which  states  that 
the  maximal  flow  value  in  a  network  is  equal  tc  the  minimal  cut 
value  is  apparent  from  the  algorithm  Just  outlined,  since  the 
at  which  the  computation  terminates  occurs  wnen  some  cut 
in  the  network  is  saturated  by  fiows  in  the  right  direction. 

We  also  point  out,  for  whn  t  it  may  be  worth ,  that  the 
assertion  in  the  footnote  on  page  9  can  now  be  verified.  To  see 
this,  let  us  define  a  3tate  of  the  system  to  be  any  flow  cor¬ 
responding  to  a  basic  partition,  where  the  arcs  of  i  are 
saturated  and  those  of  Z  have  zero  flow,  i.e.  a  state  1b  any 
flow  that  might  conceivably  be  reached  in  the  computation. 

There  are  only  a  finite  nuraoer  of  states,  Bay  Thus, 

if  we  start  with  any  state  <3^  and  randomize  among  tle3  at  each 
step,  there  is  a  positive  probability  p^  of  reaching  a  3tate 
satisfying  the  maximal ity  criterion  within  t  steps,  3ince  the 
poasibl 1 ltles  at  each  step  Include  the  one  which  would  be 

selected  in  the  perturbed  problem.  Let  p  -  min  p.  >  0.  TVien 

1  1 

the  probability  of  not  terminating  is  11m  ( 1— p)™  -  0. 

m  ->•  ao 

For  anyone  who  may  be  Interested  in  computing  maximal  flows 
in  networks  by  this  method,  we  do  not  recommend  altering  the 
capacities  as  described  in  this  section.  All  empirical  evidence, 
both  in  this  particular  problem  and  with  the  simplex  algorithm 
In  general  (of  which  this  is  a  special  case),  Indicates  that 
the  insurance  bought  In  this  way  Is  not  worth  tne  effort. 
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